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ABSTRACT 


The fractal nature of non-statistical fluctuations in the density distribution of singly charged particles produced in *”S— 
Ag/Br interactions at an incident momentum of 200 A GeV/c has been investigated under the frame work of Hwa’s 
multifractal moments, Takagi’s Multifractal moments and multifractal detrended fluctuation (MDFA) analysis. The 
experimental data have been collected by using the nuclear photographic emulsion technique. All results obtained 


experimental data is analyzed and have been compared with the simulated results. 
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1. INTRODUCTION 


As the size of the experimental data is finite there must be present fluctuation in the distribution of produced particle 
spectrum in high energy heavy ion collision (AB collision) which is statistical in nature. But along with this there also 
exists another type of fluctuation which arises due to some dynamical reasons. By taking an average over a comparatively 
large sample the statistical fluctuation may be minimized by a substantial amount. At the time of averaging, the dynamical 
components are also averaged out, as a result in the final state the distribution become smooth. There exist many statistical 
techniques by the use of which one may obtain information about the physics of the dynamical nature of the fluctuation 
present in the distribution of produced particle spectrum. In this field it is well known that the density fluctuation has self- 
similar multifractal properties which may have resulted due to some kind of scale invariant dynamics. Evaluating 
appropriate moments of the distribution and by examining how they depend on the phase-space interval size (5x) one can 
draw information about the nature of the dynamical fluctuation. According to the theory [1-3] and from the experimental 
results [4-8], it is established that the self-similarity in density fluctuations should lead to a power-law scaling behaviour of 
the 8X -dependence of multiplicity moments. Such scaling laws can further be utilized to extract universal fractal 
properties of the underlying distribution and its fluctuation. Efforts have been made to interpret the observed scale 
invariances in terms of the random cascading model, phase transition or more conventional phenomenon such as the Bose— 
Einstein correlation, but each with limited degree of success. Both the experimental and phenomenological status of the 


subject has been comprehensively reviewed in [9]. 
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The self-similarity of scaled factorial moment (SFM) or intermittency phenomenon is observed in our previous 


investigation [10]. We observed that in pseudorapidity (1) space the SFM (F, )of order g has been found to obey a scaling 


property like 


F,« 8x **, (1) 


For a region q=2 to q=6, where 18 known as intermittency index. Our previous investigation [10] established 


that the self-similarity of density fluctuations in one- dimensional particle distribution down to the experimental resolution 
involves the fractal structures. In the paper [10] we already presented a detail and systematic investigations under the frame 
work of Hwa’s multifractal moments and Takagi’s moments and measure the different fractal dimensions. Hence, in this 
paper I have presented some results on scaling behaviour of multifractal moments as functions of§X , and has been 
examined by using the different techniques [1, 11- 15]. All the techniques have been described shortly in the methodology 
section. Results obtained from these three different formalisms of multifractal analysis have been compared to the extent 
possible. **S—A/Br events at the same incident momentum and having identical multiplicity distribution as the 
experimental one have been simulated by using the computer code FRITIOF based on the Lund Monte Carlo model [16]. It 
is used to compare experimental in case of first two techniques. In order to compare the experimental results in case of 
third technique I prepared an identical set of simulated data obtained from UrQMD code [16] as experimental one and 
compare the experimental results with the simulated results. The major objectives of this paper are to establish the presence 
of multifractal characteristics of the experimental data on density fluctuations beyond those arising from statistical noise in 
terms of MFDFA moments, and compare all the results obtained from first two techniques how far they agree with each 


other as well as to extract relevant fractal parameters. 
2. EXPERIMENT AND DATA 


Nuclear photographic emulsion plates (Ilford G5) were irradiated horizontally with a beam of *’S nuclei at an incident 
momentum of 200 A GeV/c obtained from the super proton synchrotron (SPS) at CERN. After development, washing, 
mounting on glass plates, and drying, the emulsion plates were volume scanned with the help of Leitz Metalloplan 
microscopes at a total magnification of 300x. For minimize the error and biasness the scanning was done by two 
independent observers. Considering a primary interaction induced by the incoming *’S projectile the total number of track 
of secondary emission was counted and the emission angle (6 ) and the azimuthal angle () of each secondary track with 
respect to the incident projectile track the number of secondary tracks were measured by an oil immersion objectives at a 
total magnification of 1500x. According to emulsion terminology, the tracks coming out of an interaction can be classified 
into four categories, namely, shower, grey, black tracks and projectile fragments. The singly charged particles moving at 


relativistic speed(B ~ 0.7) is known as shower track. The particles produced in a high-energy interaction (mostly charged 


mesons) fall into this category, and the total number of such particles in an event is denoted by n,. The grey tracks are 


mainly due to protons knocked out from the target nucleus that directly participate in the interactions. A few percentages of 
them may also be due to slowly moving mesons. The slowly moving heavy target fragments that coming out of the target 
nucleus after the interaction has taken place is termed as black tracks. The projectile fragments are the spectator parts of the 
incoming projectile nucleus. Details of the event selection criteria, classification of tracks and other aspects of the 


experiment can be found elsewhere [17]. The number of heavy tracks in an event (n,) is equal to the sum of the number of 
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black (n,) and grey tracks (n,). A cut on the number of these heavy target fragments e.g., mn, = 8 in each event ensured 


that a subsample of **S-interactions only with Ag or Br nuclei has been considered. Here I consider only those interactions 


for which the number of spectator projectile fragments in an event (7,,) with charge Z > 2 is equal to zero, which enables 


us to choose only those interactions for which total fragmentation of the incident *’S nucleus has taken place. Following 


these criteria, only 200 central and semi-central **S—Ag/Br interactions was selected for further analysis. The average 


number of shower tracks for the considered sample of events was,<n, >= 217.79 + 6.16, and the present analysis is 


confined only to the shower tracks. 


The rapidity variable, defined as 


Is additive in nature under Lorentz transformations and is used to locate a particle in a one-dimensional phase 
space. Here E and p, are, respectively, the energy and longitudinal component of linear momentum of the particle. As 


energy and momentum of the emitted particles is very difficult to measure in case of emulsion experiment it is convenient 


to replace the rapidity variable by the pseudorapidity variable, as 
y =—Intan = 
2 


Here in comparison with the total energy the rest energy of a particle was neglected, as is the case for most of the 
charged mesons produced in high-energy interactions. There are many limitations in Nuclear emulsion experiments but it is 
superior to other experiments in one respect that here we may get a very high angular resolution (= Imrad), and this 
advantage can be exploited, where distribution of particles in small phase-space region is to be analysed. To overcome the 


problem of shape dependence of the emitted particle distribution one can replace the phase-space variable (sayY] ) with a 


cumulative variable Xn [19] defined as 


{pman 


x — Mnin 
n Tmax 


Joeman 


N min 


Here, Ninin 20d Nya are the minimum and maximum value off], and p(y)is the single-particle density 
distribution in terms of1) . Irrespective of the basic phase-space variable from which it is derived, density distribution in 


terms of the cumulative variable is always uniform in between Y =O andy =1. Though the entire analysis on multiplicity 
moments will henceforth be preformed taking Xm as the basic variable, we shall continue to call the corresponding space 


dn, 
an 


as 1] -space. The pseudorapidity density distribution p (7) ) = N,| (—*), for the shower tracks was obtained, N,, (= 200) 


being the total number of events in our sample. The distribution of p (1) is well represented by a Gaussian distribution. As 
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mentioned earlier, the experimental results have also been compared with those obtained by analysing events generated 
with the computer code FRITIOF based on Lund Monte Carlo model [21] for high-energy AB interaction in case of Hwa’s 
multifractal moments and Takagi’s multifractal moments and in case of MFDFA moments we used UrQMD code [16]for 


simulated data. 


3. METHODOLOGY AND RESULTS 
3.1 Hwa’s Multifractal Moments 


It is well known that with the help of fractal geometry it is possible to characterize the distribution and it is well observed 
in the investigation from the intermittency analysis of the present set of data [10]. According to the theory of multifractality 
the density of final state hadrons should follow a scaling property with the phase-space partition size, and the scaling 


properties should be different in different regions of phase space. 
Due to the finiteness of average shower track multiplicity <n, > multiparticle production phenomenon in high 


energy nucleus-nucleus collision suffers an acute problem. The statistical fluctuation is very large for finite <n, >of the 


frequency distribution and its moments. As the bin size gets smaller, the problem of statistical noise arising out of the 
growing presence of empty bins also requires special attention. The multifractal moment (G-moment) of order gq, also 


known as frequency moment, was introduced for an event as [11] 


M 


C= 
‘ para, 


my 


O(n, - 9) (2) 


Here M is the total number of intervals into which the entire phase space has been divided, Ni is the number of 


M 
particles in the j" bin of i event, (1, ); is the total number of particles in the i" event, [(n,), = > n, |» and © is a step 


j=l 
function for integer as well as fractional g as defined in [11]. The step function is taken care of in order to account the 
effect of empty bin. The theory of fractals demands that a self-similar behaviour in dynamical component should present in 


the density fluctuation, the G-moments should exhibit following scaling behaviour: 
G,(5x,) « 6X, 5: BY, 0. (3) 
Here 7(q) is called the mass exponent. The fractal behaviour cannot be extracted in the strict sense as the partition 
number M remains finite, as the limit 5X, — Ocannot actually be reached and. But one can still obtain significant results 


by examining the scaling properties of the G-moments in a region where dy, ‘i is of the order of the phase-space resolution 


permitted by the experiment. Taking the vertical average of G-moments over the sample of events under consideration, 


one can determine the event space average of mass exponent as 


O(< InG, >) 
o(In M ) 


<™(q) >=- (4) 
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Figure 1: Variation of G-Moments with Phase-Space Partition Number for *°S— 
Ag/Br Interaction at 200 A GeV/c: (a) Experiment, (b) FRITIOF and (c) 
Random Number. In all Diagrams the Continuous Lines are Drawn Simply by 
Joining the Points. 

In the papers [1, 2 & 11] give a systematic and brief descriptions to determine various parameters related to 
multifractal characteristics of density fluctuation. For different g Values of <InGg > have been plotted against InM in Fig 1. 
The experimental results, the FRITIOF simulated results as well as the results obtained by random number generation are 
shown in the diagram separately. From the above graph, it is very clear that for all three sets of data <InGq > linearly 


depend on InM in accordance with equation (3), increasing for g < O and decreasing for g > Ishowing a saturation effect in 


the large InM region. This saturation effect is may arise due to finiteness of <n, >. 


4.0, 
(a) eA <g> = [ (b) wf 


SOOT VI UD ODO nue 1.0} 7 
0.0F Soq> osl oy 
Fo, : 
06 Va A 


-4.0F+ / 


io 

. 

oF 0.0 ; 

-8.9 4 n n n J ns n n 1 A 4 
-6 os 0 3 q 6 06 08 10 12 14 16 18 


<O,> 

Figure 2: (a) Event Averaged Mass Exponents <z(q)> and Lipschitz— 
Holder Exponents <a, > Against q for *°S_Ag/Br interaction at 200 A 
GeV/c. Data Points Represent the Experimental Values and the Lines 
Represent the Corresponding FRITIOF Predictions. (b) Multifractal 
Spectral Function for both the Experiment and FRITIOF. The Solid 
Curve with Points Represents Experimental Results and the Dotted 
Curve Represents the FRITIOF Prediction. The Straight Line is the 
< f(Q,) >=<a, > Tine 


Following equation (4) the <7(q)> values can be determined for each g from the best linear fit of the nV 


dependence of <InGq >. The event space averaged multifractal spectral function 
< f(Q,)>=q<Q,>-<tT,>, 
Is introduced through a Legendre transform with the help of Lipschitz—Holder exponent a, , that is defined as 


O<T,> 


<Q, >= 
q oq 
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Since a derivative is involved, it is necessary to determine <a, > for small incremental changes in gq, especially in 
the neighbourhood of g = 0, where <f(a, )> has its maximum. Fig 2(a) shows the variation of <z,> and <a,> values for the 
experimental and FRITIOF data with g. Here no such significant difference is found between the experimental and 
simulated results which are contradictory to the case of intermittency analysis [10]. A smooth and stable multifractal 
spectral function <f(a, )> has been obtained both for the experiment and for the FRITIOF. Multifractal spectral function <f 
(a, )> for both experiment and simulated data are plotted against <a,> in Fig 2(b), and both satisfy the general 
characteristics mentioned in [1, 3] such as (i) <f(a, )> is a function of <a, > that is concave downwards, (ii) has a peak at 
<a >, and (iii) the straight line <f (a, )> = <a, > tangentially touches both the spectra around <a\>, because <f (a,)> = 
<a,> and <f’ (a,)> = 1. The region above the <f (a, )> = <a, > line has no physical significance and corresponds to an 
unphysical region. The wide distribution in <f (a, )> and not a delta function peaked around ap observed in the diagram 
confirms multifractal nature of the density fluctuation in each case. The region of the spectrum where the value of <a, > 
less than unity corresponds to the dense region of the spectrum and for the region for which the value of <a, > is greater 
than unity is corresponds the sparse region of density distribution. Both experimental and simulated maximum values of <f 
(aq )> are very close to unity, indicating that the empty bin effect particularly in the higher resolution region is marginal in 
the present case. The FRITIOF simulated spectrum is wider than the experimental one. We already know for a pp 
interaction a similar feature has been observed while UA1 data were being compared with GENCL and PYTHIA 
predictions [4]. In the case of AB interactions multifractal characteristics were observed both in the experimental data as 


well as in the Monte Carlo predictions based on a simple stochastic model [8]. 
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Figure 3: (a) Experimental Values of the Intermittency Indices (9, ), (q — 1 — <t(q)>) and (q -1- 
<z(q)>*") are Plotted Against g for *S-Ag/Br Interaction at 200 A GeV/c. The Straight Lines are 
best Linear fit to data. (b) Experimental Values of d,/d, Obtained both from SFM (Solid 
Diamonds) and G-Moments (Open Diamonds) are Plotted Against g. The Continuous (Solid and 
dotted) Lines Represent the Corresponding Levy Law Prediction using Equation (9). 


For the finiteness of the average number of charged particle, the G, moments contain statistical contribution (G? ); 


this statistical part can be determined by distributing ns particles of an event randomly within0 < Xn <1. It should be 


remember that in this process short-range correlation among the particles, if any, is destroyed. The dynamical contribution 


(Go") can then be extracted after eliminating the statistical one. In [11] it has been shown that, for a trivial dynamics the 


dynamical part of <r, >, denoted by <7(q)>“", should be equal to (q — 1). The deviation in <7(q)>°" from (q -1) is from the 


nontrivial dynamical contribution. As all three G-moments G,, G* and Ge obey their respective power laws, the 
following relation may be obtained: 


dyn. 


<1(q)>"=<(q)>-<t(qg >" +9q-1 (5) 
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whereas, the intermittency index g, introduced in equation (1), can also be connected to <7(. qo as [11] 


<t(q) >" -q+1~9,. (6) 


In Figure 3(a) the variation of g, values along with the (g — 1 — <r(qg)>) and (q- 1 - <xq)>" 


) values against q is 
shown. From the diagram it is clear that the g, values differ from the respective (q — 1 — <1(q)>") values by very small 
amount. The difference in their values may probably be attributed to the different ways of defining SFM and G-moment. 


The generalized Renyi dimensions, denoted by D, , may be obtained from the intermittency indices as 
p, =1-2_ (7) 


Therefore, with the help of equation (6) one can demand 


dyn 


D w tg | (8) 
q q-1 


Simultaneously the anomalous dimensions are defined as 
d,=D-D,,. 


Where, D is the topological dimension of the supporting space. D = |stands for the one dimensional analysis. The 
Levy index ( ) is a very useful parameter to classify the properties of universal multifractals. For multifractal nature it is 
found that the value of is lying between O< < 2. Levy index indicates the degree of multifractality as well as estimates 
the cascading rate in self-similar branching process [16]. The Levy index (_ ) can also be utilized to highlight the possible 
mechanism of particle production. Such a characterization of multifractality is possible if the underlying density 
distribution can be described by a Levy stable law. Under a Levy law approximation, using anomalous dimensions one can 


determine the value of from the following relation [21]: 


Ps 
d, 2° =2 @-1 


In Figure 3 b. Experimental values of d,/d, obtained both from SFM (solid diamonds) and G-moments (open 
diamonds) are plotted against g. The lines correspond, respectively, to = 1.15 for the exact values of g, and to = 1.25 
for the approximate values. If were equal to 2, the Levy distribution would have transformed into a Gaussian one. 
Minimum fluctuation in the self-similar branching processes is expected under such condition. On the other hand, the value 
of dq/d2 become independent of order for = 0. This condition corresponds to mono-fractals and maximum fluctuation. 
Under such condition one may expect a second-order phase transition. In our case neither of the above two conditions is 
satisfied by the values obtained. One can see that the present values obtained from two different sets of parameters, are 
very close to each other. The fact that >1.0, demands the presence of wild type of singularities arising out of fluctuations 
in the density distribution which is non-Poisson like in nature. Here, it can be concluded as the mechanism of particle 
production in the present case of **S—Ag/Br interactions may be described in terms of a possible non-thermal phase 
transition during the cascading process. On the other hand, a value of 0 < = < 1.0 would have indicated soft bound 
singularities, that can be related to a thermal phase transition interspersed in the cascading process. It should however be 


mentioned that the present values of Levy index are less than a previously obtained value ( = 1.6) based on a set of 
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combined data on AB, pA, e+e— and _p interactions [21], but are well within the limit allowed by the Levy law description, 
and do not necessarily warrant a thermal phase transition to occur during particle production. 

3.2 Takagi’s Multifractal Moments 


In order to overcome the difficulties of finiteness of the number of charged particle multiplicity an alternative approach has 


suggested by Takagi [12]. According to his methodology a new set of multiplicity moments for g > 0 are suggested as 


7, (8x4) =In > (p,)" (10) 


i=l j=l 


T,, is known as Takagi’s multifractal moments and are not affected by the finiteness of n,. Here, pij (= nij /K) is the 
normalized density function, K is the total number of particles produced in Nev interactions and nij is the same as equation 
(2). Takagi’s method is based on two assumptions, (i) over the considered phase-space interval the density function p is 


uniform, and (ii) the multiplicity distribution Pn does not depend on the location of the interval $y & 2 The above two 


conditions are found to be valid in the present case where X7 has been used as a phase space variable. According to the 


theory of multifractals, Tq( dy ; ) should be a linear function of the logarithm of the resolution R( $y ‘ ) as, 
T,6X,,) =A, + B, In ROy,,): (11) 


Where Ag and Bg are two constants which are independent of q. If the linear dependence like that of equation (11) 


is observed over a large range of R(dy _ ), following Takagi’s method the generalized dimensions can once again be 


calculated as - 
D, =B,(q-1) (12) 


For g = 1, one may either take the appropriate limit [22], or may introduced the entropy function defined as 


Ney 
s(x,)=->'> p,lnp, (13) 


i=l j=l 


The entropy function is also defined as, 
S(8x,,) = —D, In R(6X,,) + const.» (14) 


Where, D1 is termed as information dimension. If the number of events Nev is quite large one will get 


yy (p,)" =<n,>K,,<n>) - (15) 


i=l j=l 
Where average bin multiplicity <n> = K/(M. Nev), and therefore, it can be write 


Inn, = A, +[(q¢-DD, +1ln@y,) (16) 


For the simplest choice of R( $y ‘ ), = 8x a Replacing Sy ‘ with <n>, the generalized dimensions can now be 


obtained following the relations: 
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Inn! =A,+[(q-DD, +1]In<n> (17) 
for g 22. For g=1, 
<ninn>/<n>=C,+D In<n>. (18) 


Proceeding in the same way as [12], for a symmetric interval about the central value ) 0.5 of the distribution, 


The values of <n In n> and In<nq > are calculated with increasing width of the interval about the central value ( 0.5) 


of the distribution according to the procedure mentioned in [12]. Our results on Takagi’s method of multifractal analysis 
have been graphically shown by plotting <n In n>/<n> against In<n>, in Fig 4(a)-(c), respectively for the experimental 
values in both 7 and ¢-space and the FRITIOF prediction in 7-space. Similarly the variation of In<n‘> against In<n> is 
shown in the Fig 4(d)-(f). From the slopes of best linear fit to data values of the plot <n In n>/<n> vs In<n> the 


information dimension (D;) have been calculated respectively as D, = 0.973 + 0.0014 ,0.972 + 0.002 and 0.979 40.002. 


Values of generalized dimensions for g > 2 have been obtained 


60 r 
(b) Experiment 


phi space 


(a) Experiment 
cta space 


Bn 
° 
' 


<n In n=/<1> 


[ (c) FRImIOF 
i 
F 


2.0 s* a ck 
0.0 4 1 4 4 4 4 4 n 

20 30 40 5.0 60 370 30 40 4O 60 25 3.0 40 50 5.0 
15.0 = r : | r q=7 
Ost ot q=7 (@) Experiment re q (f) FRITIOF en : 
eta-space a=6 pnispace q- eta-space i q=6 
244 Pa A + u 
30.0 ~* ar oot es q-§ ? -5 
- oe G=3 _ gee 4 ae > ee 7 4 

Vv - a - 
ae See | Seccememe Ff a=? 
Oe" q-2 94 4 OH q-2 =2 
0.0 1 1 1 4 ' ' 1 ' ' ' ' 
2.0 3 40 5.U 6.0 20 an 40 50 60 2o 30 40 50 59 

Incn> Ini <tt> Inen> 


Figure 4: Plot of Takagi’s Multifractal Moments for **S—Ag/Br Interaction at 200 A 
GeV/c. In all Diagrams the Straight Lines Represent best Linear fit to data 


From the best linear fit of In<ng> values against In<n> as shown in Fig 4(d)-(f). For comparison, Dg values of 


different orders obtained from Takagi’s generalized moments are plotted against g >2in Fig 5, together with those 


obtained from the intermittency indices (g, ) and from the dynamical part of Hwa’s multifractal mass exponents <1(q)> ayn; 
respectively, making use of equations (7) and (8). With increasing g in general we find a monotonous deceasing trend in 
the Dq values. However, the Dg values from Takagi’s method exhibit steepest fall, whereas those obtained from the 
intermittency indices decrease at the slowest rate. Probably because of the different ways of defining the multifractal T- 
moments, and due to the reason that in Takagi’s method no attempt has been made to separate the nonstatistical 
contribution from the statistical one, in this case at large gq the Dq values differ significantly from those obtained from the 


SFMs and G-moments. However, for a simple Poissonian multiplicity distribution within a given interval $y Pe the D, 


values would all have been equal to unity. Any deviation in their values from 1.0 would thus provide us with a measure of 
nonstatistical fluctuation. This has been found in all the methods described above for characterizing multifractality in 


density fluctuations. 
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Figure 5: Experimental Values of the Generalized 
Dimensions D, for *S—Ag/Br Interaction at 200 A 
GeV/c. Solid Lines are drawn to Guide the Eye. 
On the basis of the fact that only Bernoulli type of fluctuations are responsible for a transition from monofractality 


to multifractality, Bershadski [23] gave a thermodynamic interpretation of the observed results in terms of a constant 


specific heat c: 


pap (19) 
q-1 


A monofractal to multifractal phase transition corresponds to a gap in the value of c from c = 0 to a nonzero finite 
multifractal specific heat. By plotting Dg against In gq — 1) we can, therefore, obtain the value of specific heat from the 
slope of the best linear fit. Such a plot can be found in Fig 6, both for the experiment and FRITIOF simulated values. In 


X, “Space a strict linearity is not seen over the entire range of g under consideration, and a linear fit in the range g = 2-5 


resulted into c = 0.329 + 0.061 =1/3. On the other hand, 
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Figure 6: Graph to Determine the Multifractal 
Specific Heat for *S—Ag/Br Interaction at 200 A 
GeV/c, from the D, Values Obtained by using 
Takagi’s Method. The Experimental data Points in 
n-Space are Simply Connected by a Continuous 
Line, whereas, for both the FRITIOF Prediction in 
y-Space and the Experimental data Points in ¢- 
Space, the Straight Lines Represent best Linear fit. 


best linear fit of FRITIOF data in X “Space over q = 2-7 results in c = 0.066 + 0.009, which is much smaller than 
the corresponding experimental value. On the other hand, in y, 9 “Space linear fit once again over the entire range of g (= 2- 


7) gave us a much smaller value of c (= 0.086 + 0.008). The c-value in y, » “Space is significantly smaller than a previously 
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obtained value (c = 1/4) based on an analysis in the azimuthal angle space of similar **S-ion induced experiment at same 
incident energy [24]. As a probable reason for this discrepancy it may be pointed out that the analysis presented in [24] has 
been performed over a set of events with minimum bias, that possesses a much wider range of impact parameter values, 
and also has a much smaller value of average shower multiplicity as well as smaller statistics than the present set of 
experimental data. Here the value of c in ¢-space is also smaller than the predictions of another analysis (c = 0.56) on 


simulated high multiplicity single jet events in the azimuthal angle space [25]. 
3.3 MFDFA Moments 


In many cases the local singularities of particle density functions p (5x )can be described by a power law like, 
jim p (8X ) oc (8x )’* (20) 


and may be interpreted as a self-similar (or self-affine in two dimension) fractal object with dimension Dy . Now 
each phase-space interval will have its own 6X dependence. If all of them are characterized by a same Dr, the distribution 
is considered as a geometrically monofractal object. If the exponent varies at different X positions then it is a multifractal 
one. The multifractal detrended fluctuation analysis (MF-DFA) has so far not been used extensively for multiparticle 
emission data analysis [15]. The detrended fluctuation analysis (DFA) method [13] is a very useful technique for the 
determination of (mono) fractal scaling properties and the detection of long-range correlations in noisy and stationary time 
series data [14]. Kantelhardt et al. [15] have advanced the DFA method for nonstationary and multifractal series, and the 
generalized DFA is said to be the multifractal DFA (MF-DFA) method. The DFA and the MF-DFA methods are very 
standard techniques for the time series data analysis, for completeness we provide a brief description of the methods in the 
following section. Let x: k = 1, 2,..., N be a fluctuating series (signal) of length N. The DFA/MF-DFA procedure 


consists of the following five steps: 


Step-1: Determine the Profile 
V(i)=> [eg - (x) 121.2000 (21) 
K=1 


Where (=) = >) Hs is the mean value of the analyzed signals. 
K=1 


Step 2: Determine the Profile 


Divide the profile Y(i) into N, =int(N/s) non-overlapping segments of equal lengths. One has to choose the § value 


depending upon the signal length. In case, the length N is not a multiple of the considered scale parameter s, the same 
dividing procedure is repeated starting from the opposite end of the signal. Hence, in order not to disregard any part of the 


signal series, usually altogether 2N, segments of equal length are obtained. 


Step 3: Calculate the local trend for each of the 2N, segments. This is done by least-square fits of the segments (or 


subseries). Linear, quadratic, cubic or even higher order polynomial may be used to detrend the signal, and accordingly the 
procedure is said to be the MF-DFA1, MF-DFA2, MF-DFA3, . . . analysis. Let, y, be the best fitted polynomial to an 


arbitrary segment V of the signal. Then determine the variance as 
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F,s)=>—> (IW -)s+il-y,@P (22) 


s i=l 


For N = /, 2,...... N, and V =N,+I, ....... 2N, is given as, 


F3(,8)=—5 (YIN W -N,)st+i]—y,@Y (23) 


i=l 
Step 4: The DFA function F is defined by averaging F'*(V , s) over all the 2N , Segments, i.e. 


a 


2N, 
at DOE al a (24) 


s isl 


F(s)= 


On the other hand, the q"" order MF-DFA function F. 5) is defined as 


2N l/q 


F(s)= 3D {F?W,8)) 05) 


s isl 
for all g #0 and for g = Othe above definition is modified to the following form: 


F(s) = exp [ar Eire.o| (26) 


s i=l 

In general, the order parameter g can take any real value. Note that for g = 2 the MF-DFA function reduces to the 
standard DFA function as defined in Eqn. (23). 

Step 5: Now varying the scale parameter (S ), one can study the scaling behaviour of the detrended fluctuation 
functions. If the signal series x, possesses long-range (power-law) correlation, F(s) as well as F(s) for large values of 


S would follow a power-law type of scaling relation such as 
F(s)=s" and F,(s) « s* (27) 
The exponent H is nothing but the well known Hurst exponent that can be related to the fractal dimension Das: 


D=2-H. According to the value of H, x, is considered as long-range anti-correlated if 0 < H < 0.5; uncorrelated if H = 


0.5 and long-range correlated if H >0.5. On the other hand, the h(q), is termed as the generalized Hurst exponent [15]. 


For stationery time series h(q = 2) = H i.e. the well known Hurst exponent. For a monofractal signal h(q) is independent 
of g since the variance F’(v,s) is identical for all the subsignals, and hence Eq. (25) yields identical values for all q. 
Here the function F(s) can be defined only forsm+2, where ™ is the order of the detrending polynomial. Moreover, 
F(s) is statistically unstable for very large S(N/4). If small and large fluctuations scale differently, there will be a 
significant dependence of h(q)onq. Whereas for positive values of ¢, F,(s) will be dominated by the large variance 


which corresponds to the large deviations from the detrending polynomial, for negative values of q major contributions of 


F(s) arise from small fluctuations from the detrending polynomial. Thus for positive (negative) values of g, F(s) 
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describes the scaling behaviour of the segment with large (small) fluctuations. Knowing the values of h(q) one can easily 
estimated the multifractal mass exponentt (q) and the multifractal singularity spectrum f(Q). It is very well known that 


the h(q) is related to T(q) through 
t(q) = qh(q)-1 (28) 
A non linear T(q) spectrum demands the existence of multifractal nature of the analyzed time series data. For a 


monofractal system T(q)should be a linear function of g. The singularity strength function q@(qg)and the multifractal 


spectrum f(q) are connected via Legendre transformation [26, 27]: &(qg) = OT ; / 0q as, 


f(Q) =qa -T(q) (29) 


The generalized fractal dimension D , may be derived from the mass exponent T (q) as: 


D. me (30) 
q q-1 


Here these methods are employed to the single-event1| distributions of produced charged mesons in high-energy 
nucleus-nucleus collisions. The scale parameter S is nothing but the number of partitions in T] -space. In Fig. 7, I plot the 
averaged DFA function < F'(s) > with phase space partition number (scale) S for the *°S - Ag/Br data. From Fig.1 it is 
clear that < F(s) >for large values of S would follow a power-law type of scaling relation. In this analysis I vary S from 


N to N/4, where N is the multiplicity cut. In the region 8 < 5 < 24 the scaling relation (27) holds good. As the scaling 


behaviour holds well for a fairly large value of S in all cases, one 


In<F(s)> 
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Figure 7: The Variation of Event Averaged 


DFA Fluctuation Functions < F(s)> with 
Scale Parameters for the *°S - Ag/Br 
Interaction at 200A GeV/c. The Solid (dotted) 
lines Represent the best Fitted Straight Lines 
to the Experimental (UrQMD) data. 


May expect that the signal series x, should possess some kind of long-range (power-law) correlation for both 


experimental and simulated (UrQMD) data. Here I find that H = 1.474 + 0.017 (experiment) and H =1.585+0.017 
(UrQMD). In all the cases the value of H greater than unity again establish the fact that there must present some long range 


correlation in both experimental and simulated (UrQMD) data. Using the relation Dr = 2 — H, I estimate the values of 
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fractal dimension (Dy) as Dr = 0.526 + 0.017 (experiment) and Dr = 0.415 + 0.017 (UrQMD). It is well known that, for a 
fractal object the fractal dimension D, is less than Dy, the topological dimension of the supporting space. For one 
dimensional analysis the value of D; = 1. The deviation of D; from Dy is a measure of the degree of fractality. All the 


values of Dr obtained from this analysis are much smaller than unity and hence my DFA results demand that the 1 


distributions in the interactions are highly fractal. It is also clear that the UrQMD model possesses fractal dynamics that is 
apparently identical to those of the experiments. In our previous multifractal analysis using the data of *S -Ag/Br 
interaction we observed similar behaviours. All these observations (present and previous) suggest that (multi)fractality gets 
weakened with increasing multiplicity. In Fig. 8, the variation of event averaged MFDFA fluctuation functions with scale 


parameter s is shown for both experimental and simulated data. 


*S_Aq/Br interaction at 200A GeV/c 


In<Faq(s)> 
In<Faq(s)> 3 als) 
Experiment 2h 
Or 
a} 
1.0 1.5 20 25 3.0 35 4.0 1.0 40 
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Figure 8: The Variation of Event Averaged MFDFA Fluctuation Functions with Scale 
Parameter s. Left Panel: Results obtained from Experimental data and Right Panel: 
Results Obtained from UrQMD data. The Lines Joining Points are shown to Guide the 
Eye. 
From the variation of < F q > against the scale parameter Sit is obvious that there are no saturation effects at 


higher value of S like G, moments [10] and hence one may conclude that < i > function presented here are not 


significantly influenced by the finiteness of the event multiplicity. In Fig. 9 the experimental and simulated values of Hurst 


exponent /(q) for *°§ - Ag/Br interaction are plotted against order g. The values of h(q) is determined from the linear 
portion4 < § <12 of <F, >vs S graph. Forg S$ —lregion the values of h(q) decreases slowly and there is a sharp 


fall in the g =-1 to q =1 region and then h(q) tends to saturate in the region g = 1. 
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Figure 9: Left Panel: The Variation of Event Averaged , Functions with scale 
Parameter s for *’S - Ag/Br Interaction at 200AGeV/c. The Lines Joining Points are 
Shown to guide the Eye. Right Panel: The Variation of Event Averaged Mass 


Exponent *(4) with Scale Parameters for °°S - Ag/Br Interaction 


In Fig. 9. I have also graphically represented how the values of t(q) for the experimental and UrQMD data vary 
with g . Unlike in the case of intermittency analysis, here no significant difference between the experimental and simulated 
results is observed. Mass exponents t (q) were derived from the fluctuation functions for q values between —5 and +5 and 
plotted against the gq values. A nonlinear <7 (q) > function means multiple scaling, which requires a hierarchy of scaling 
exponents (multiscaling) in order to accurately represent the scaling property. The degree of non-linearity of <t(q) > 


function can give an idea about the degree of multifractality. These nonlinear functions have convex downward facing 
plots, with the degree of convexity reflecting the level of heterogeneity in scaling exponents. A smooth and stable 


multifractal spectral function f(@)has been obtained both for the experiment and for the UrQMD for *°§ induced 


interactions. Multifractal spectra are plotted against Ol in Fig. 10, and satisfy the general characteristics [1, 2] such as (i) 


f(a)is a function of that is concave downwards, and (ii) has a peak ata, - But unlike the Hwa’s multifractal analysis 


the moments violates the condition (iii) i.e. straight line f(~) =a tangentially touches both the spectra. 


F(a) 


o4 
—=— Expenment 


—e— vramco 


Figure 10: Left Panel: Multifractal Spectrum Function F@) for both Experiment 
and UrQMD. The Straight line Represents F(A) =o 4 Right Panel: 


Generalized Fractal Dimension Dy as a Function of @ for both Experiment and 
UrQMD in Case of **S - Ag/Br Interaction. 
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The region above the line f(a) =a corresponds to an unphysical region. The fact that a wide distribution in 


f(Q) and not a delta function peaked around Ol, has not been obtained. In the present analysis the peak values of the 


spectra are shifted to a higher value of Ol. All observations confirm the multifractal nature of the density fluctuation in 
each case. The left and right side of the spectrum corresponds, respectively, to the dense and sparse regions of density 


distribution. Both simulated and experimental maximum values of f(q) are very close to unity, indicating that the empty 


bin effect particularly in higher resolution region, is marginal in the present case. The generalized fractal dimension D, 
plotted against the order q in Fig. 10, the values of generalized fractal dimensions D, determined using MFDFA methods 


are as follows: 


Table 1 
. Experiment UrQMD 
Interactions | 4 D, | 8D, D, | SD, 
2 | 0.437 | 0.0376 | 0.361 | 0.0440 
2 3 | 0.661 | 0.0178 | 0.601 | 0.0213 
S- Ag/Br 410.752 | 0.0114 | 0.698 | 0.0139 
5 | 0.805 | 0.0083 | 0.753 | 0.0102 


From the tabulated values of D, obtained from the MFDFA method it clear that all values are much smaller than 


those obtained from the previous conventional methods. It is also clear from the tabulated values that there exists 
insignificant difference between the experimental values and the corresponding UrQMD predictions in case of '°O - Ag/Br 
interactions, but there exists some discrepancies in case of *S - Ag/Br Interactions. In this present investigation a distinct 


and systematic increase in is found with the order number q. This behaviour is quite different from the behaviour in the D, 
values obtained from the conventional methods. It should mention that the D, values obtained from Takagi’s method 


exhibit a sharp fall. This is probably because of the fact that Takagi moments are not free from the statistical noise [12]. 
4. DISCUSSION 


In this analysis I presented a systematic analysis of the pseudorapidity fluctuation of charged mesons produced in *’S - 
Ag/Br interaction at 200A GeV in terms of the Hwa’s multifractal moments, Takagi’s moments and DFA and MF-DFA 
methods. To understand the underlying mechanism(s) of particle production in these interactions, the experimental results 
are simulated by using the Fritiof model in case of first two moments and UrQMD model for the third technique. From the 
above analysis it is established that experimental data along with all simulated data behave like (multi)fractal system. The 


nature of h(q).t(q)and the singularity spectra also confirm the presence of multifractality in the data as well as in the 


simulation. The nature of these spectra and the estimated values of the Hurst exponent H demand that the origin of 
fractality is of the two, three or higher order particle correlation in all cases. The MF-DFA prediction of the generalized 
fractal dimensions are consistently lower than that obtained from previous techniques like scale factorial moment analysis, 
Hwa’s multifractal analysis and Takagi’s method. In case of MFDFA moments it should noted that within the error 
margins the experimental results cannot be discriminated from their UrQMD values. The observations signify that the 
MFDFA technique like Hwa’s multifractal moments and Takagi’s moment is probably not sufficiently sensitive to the 
nature of fluctuation present in the data. A reliable method of filtering out the statistical noise from the MF-DFA function 


is needed to make the technique more effective. 
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